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Abstract 


A parametric study is presented of the buckling behavior of infinitely 
long , symmetrically laminated anisotropic plates subjected to combined 
loadings. The loading conditions considered are axial tension and com- 
pression ■, transverse tension and compression , and shear. Results ob- 
tained using a special-purpose analysis, well-suited for parametric stud- 
ies , are presented for clamped and simply supported plates . Moreover , 
results are presented for some common laminate constructions, and 
generic buckling design charts are presented for a wide range of pa- 
rameters. The generic design charts are presented in terms of useful 
nondimensional parameters, and the dependence of the nondimensional 
parameters on laminate fiber orientation, stacking sequence, and mate- 
rial properties is discussed. 

An important finding of the study is that the effects of anisotropy 
are much more pronounced in shear-loaded plates than in compression- 
loaded plates. In addition, the effects of anisotropy on plates subjected 
to combined loadings are generally manifested as a phase shift of self- 
similar buckling interaction curves. A practical application of this phase 
shift is that the buckling resistance of long plates can be improved 
by applying a shear loading with a specific orientation. In all cases 
considered in the study, the buckling coefficients of infinitely long plates 
are found to be independent of the bending stiffness ratio (D\\f -C>22 ) . 


Introduction 

Buckling behavior of laminated plates is a topic 
of fundamental importance in the design of aerospace 
vehicle structures. Often the sizing of many sub- 
components of these vehicles is determined by stabil- 
ity constraints in addition to other constraints such 
as stiffness and strength. One typical subcomponent 
of great practical importance in structural design is 
the long rectangular plate. These plates routinely 
appear as subcomponents of stiffened panels used for 
wings, ribs, and spars. In addition, long plates com- 
monly appear as subcomponents of semimonocoque 
shells used for fuselage structure, liquid-propellant 
booster tankage, and booster interstage structure. 
Buckling results for infinitely long plates are impor- 
tant because they provide a lower bound on the be- 
havior of finite- length rectangular plates, and they 
provide information that is useful in explaining the 
behavior of finite-length plates. 

Understanding the buckling behavior of symmet- 
rically laminated plates is an important part of 
the search for ways to exploit plate orthotropy and 
anisotropy to reduce structural weight of aircraft and 
launch vehicles. Symmetrically laminated plates can 
exhibit anisotropy in the form of material-induced 
coupling between pure bending and twisting defor- 
mations. This coupling generally yields buckling 
modes that are skewed in appearance, as depicted 


in figure 1. The effects of anisotropy on the buck- 
ling behavior of these compression- loaded plates is 
well-known (e.g., refs. 1 and 2). However, the effects 
of anisotropy on symmetrically laminated plates sub- 
jected to shear loadings and various other combined 
loadings are generally not as well- understood. 

Symmetrically laminated plates can have many 
different constructions because of the wide variety of 
material systems, fiber orientations, and stacking se- 
quences that can be used to construct a laminate. 
Thus, it is extremely useful for the designer to have 
a means of assessing the performance of various lam- 
inated plates on a common and unbiased basis. The 
use of nondimensional parameters can provide this 
means for assessing laminate performance. More- 
over, the use of nondimensional parameters permits 
buckling results to be presented in a concise manner 
on one or more plots as a series of generic curves that 
span the complete range of plate dimensions, loading 
combinations, boundary conditions, laminate con- 
struction, and material properties. An added ben- 
efit is that the generic curves also furnish the de- 
signer with an overall indication of the sensitivity 
of the structural response to changes in the design 
parameters. Examples of generic design charts that 
use nondimensional parameters are presented in ref- 
erences 3 and 4 for buckling and postbuckling of or- 
thotropic plates. 


The major objectives of the present paper are 
to indicate the effects of plate bending orthotropy 
and anisotropy on the buckling of plates subjected to 
combined loads, and to present generic design charts 
that indicate the buckling behavior of long plates for 
a wide range of parameters in a concise manner. The 
loading conditions considered in this paper are ax- 
ial tension and compression, transverse tension and 
compression, and shear loadings. Results are pre- 
sented herein for plates with the two opposite long 
edges clamped or simply supported. A substantial 
number of generic buckling curves that are applica- 
ble to a wide range of laminate constructions are pre- 
sented in the present paper. 

Symbols 

A n , B n displacement amplitudes 

(see eq. (8)), in. 

b plate width (see fig. 1), in. 

£>ll , D\ 2 , T>22> A36 orthotropic plate bending 
stiffnesses, in-lb 

£>26 anisotropic plate bending 

stiffnesses, in-lb 

E\, E2, G\2 lamina moduli, psi 

K x , K y . K s nondimensional buckling 

coefficients associated 
with axial compression, 
transverse compression, and 
shear loadings (see eqs. (5), 
(6), and (7), respectively) 

ra laminate stacking-sequence 

number 

N x , N y , N xy membrane prebuckling 

stress resultants (see fig. 1), 
lb/in. 

Af£ r , Ny r , N£!y values of membrane stress 

resultants at buckling, lb/in. 

p, p cr nondimensional loading pa- 

rameter and correspond- 
ing value at buckling, 
respectively 

wj\i out-of-plane displace- 

ment field at buckling (see 
eq. (8)), in. 

x, y plate coordinate system (see 

fig. 1), in. 


a, /?, 7, 6 nondimensional param- 

eters defined by equa- 
tions (1), (2), (3), and (4), 
respectively 

p = y/b, £ = x/\ nondimensional plate 

coordinates 

6 fiber orientation angle (see 

fig. 1), deg 

E 

A half-wavelength of buckling 

mode (see fig. 1), in. 

I 

v\2 lamina major Poisson’s ratio l 

4> n , 4> n kinematically admissible 

basis functions (sec eq. (8)) j 

i 

Approach 1 

The objectives of the present study are achieved ] 
herein in two different ways. In the first way, the cf- 1 
fects of plate orthotropy and anisotropy arc presented 
in an implicit manner for symmetrically laminated 
angle-ply plates. This family of laminates is denoted | 
herein by the symbols [(±0) m ] s , where 6 is the fiber 
orientation angle shown in figure 1 (given in units of | 

degrees) and m is the stacking- sequence number that | 

indicates the number of times that the plies in paren- - 

theses are repeated. This family of laminates was 1 

chosen because varying the laminate parameters 6 = 

and m encompasses a wide range of plate bending or- : 

thotropy and bending anisotropy. The term implicit f 

is used above to reflect the fact that for the angle-ply j 

laminates, as well as for all other laminates, the plate I 

orthotropy and anisotropy are implicit functions of | 

laminate fiber orientation and stacking sequence and 
do not typically vary independently. This implicit 
way of assessing plate behavior is generally the ap- = 
proach that is the most familiar to aerospace vehicle 
designers. 

A second way to characterize the effects of plate 
orthotropy and anisotropy uses the nondimensional 
parameters presented in reference 1. These non- 
dimensicnal parameters are used in the present paper 
to show the effects of plate orthotropy and anisotropy 
on buckling behavior in a somewhat more explicit 
manner. For example, buckling results are presented j_ 
in the present paper as a function of two parameters j_ 
that characterize plate orthotropy and two param- 
eters that characterize plate anisotropy. By vary- 
ing each parameter in a systematic manner, the ef- 
fects of orthotropy and anisotropy on buckling can be 
assessed independently. An important step in this as- 
sessment is to establish an understanding of how lam- 
inate fiber orientation, stacking sequence, and mate- 
rial properties affect the values of the nondimensional ! 
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parameters. Thus, results are also presented herein 
that indicate the effects of laminate construction on 
the nondimensional parameters for a few selected 
laminates. Material properties considered herein 
are representative of high-strength graphite-epoxy 
material, ultrahigh-modulus graphite-epoxy mate- 
rial, S- glass 1 -epoxy material, Kevlar 2 -epoxy mate- 
rial, boron-epoxy material, and boron- aluminum 
material. 

Analysis Description 

Often in preparing generic design charts for buck- 
ling of a single plate element, a special-purpose analy- 
sis is preferred over a general-purpose analysis code, 
such as a finite element code, because of the cost 
and effort involved in generating a large number of 
results with a general-purpose analysis code. The 
results presented in the present paper have been ob- 
tained using such a special-purpose analysis. A brief 
description of the analysis is presented subsequently. 

The buckling analysis used in the present study 
is based on the classical Rayleigh-Ritz variational 
method and is derived explicitly in terms of the 
nondimensional parameters defined in reference 1. 
Deriving the analysis in this manner inherently 
makes the resulting computer code well-suited for 
parametric studies. The nondimensional parameters 
used are given by 



DjU + 2D66 

V D\ \ D22 
Dig 

(^22) 1/4 
£>26 

(Dn ^ 2 ) 1/4 


( 1 ) 

(2) 

( 3 ) 

( 4 ) 


where A is the half- wavelength of the buckling mode 
and 6 is the plate width shown in figure 1. The sub- 
scripted D terms appearing in the equations are the 
plate bending stiffnesses of classical laminated plate 
theory. The parameters a and /? characterize plate 
bending orthotropy, and the parameters 7 and 6 char- 
acterize plate bending anisotropy. The parameters 
defined by equations (2) through (4) depend only on 
the plate bending stiffnesses, whereas the parame- 
ter a depends on the buckle aspect ratio A/6 also. 


^S-glass: Trademark of the Ferro Corporation. 

^Kevlar: Trademark of E. I. duPont de Nemours & Co., Inc. 


Without loss of generality, and as a matter of conve- 
nience, the nondimensional parameter (D\\/ D22) 1 / 4 
is used in the present study in the place of the pa- 
rameter a to assess plate bending orthotropy. 

In addition to {D\\j D22) 1 ^ , P, 7? an ^ ^ three ad- 
ditional nondimensional quantities are used to char- 
acterize buckling resistance, and they are given by 


K x = 

N?b 2 

( 5 ) 

ir 2 \J D\\D22 

Ky = 

Ny r b 2 

( 6 ) 

TX 2 D 22 

Ks = 

Nxyb 1 

( 7 ) 

^ {d u d% 2 ) 1/4 


The quantities K x ,K y , and K s are referred to herein 
as the axial compression buckling coefficient, the 
transverse compression buckling coefficient, and the 
shear buckling coefficient, respectively. Each of the 
loading conditions associated with these buckling co- 
efficients is shown in figure 1. The positive- valued 
compressive stress resultants and the shear stress re- 
sultant corresponding to the loadings are denoted by 
N Xi Ny , and N xy , respectively. The stress resultants 
appearing in the buckling coefficients correspond to 
critical values associated with the onset of buckling. 

In the buckling analysis, the infinitely long plates 
are assumed to have uniform thickness and material 
properties that do not vary along the plate length and 
width. In addition, the uniform biaxial and shear 
loadings and the boundary conditions do not vary 
along the plate length. Under these conditions, in- 
finitely long plates have periodic buckling modes that 
exhibit either inversion symmetry or inversion anti- 
symmetry with respect to a given reference point. 
The buckling mode depicted in figure 1 corresponds 
to either inversion symmetry or antisymmetry de- 
pending on the selection of the reference point. The 
presence of the periodicity and inversion symmetry 
or antisymmetry in the plate buckling mode aids in 
simplifying the buckling analysis. 

An important characteristic of an infinitely long 
plate that possesses a periodic buckling mode is that 
a basic repetitive buckle pattern can be identified 
from which the overall buckling mode can be con- 
structed by axial translation of the pattern. The se- 
lection of the basic repetitive buckle pattern is not 
unique, and thus one may choose it to simplify the 
buckling analysis further. In the present study, the 
basic repetitive buckle pattern was chosen because 
it possesses inversion symmetry about its geometric 
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center. This simplification is important in that it 
permits the elimination of superfluous terms in the 
kinematically admissible series used to describe the 
buckling mode, which in turn simplifies the ensuing 
analysis and reduces the number of computations re- 
quired to obtain an accurate result. The mathemat- 
ical expression used in the variational analysis to de- 
scribe the buckle pattern is given by 

N N 

W N (£, n) ~ ^ A n $»(/?) + cos 7r£ (8) 

ri=l n— 1 

where £ = xf A and rj = y/b are nondimensional 
coordinates (see fig. 1), wjy is the out-of-planc dis- 
placement field, and A n and B n are the unknown 
displacement amplitudes. In accordance with the 
Raleigh-Ritz method, the basis functions $ n (v) and 
^n(v) are required to satisfy the kinematic boundary 
conditions on the plate edges at 7] = 0 and 1. For the 
simply supported plates, the basis functions used in 
the analysis arc given by 

(ji) = sin 2n7n7 (9a) 

^ n(v ) = sin(2n - (9b) 

for values of n = 1, 2, 3, . . . , N. Similarly, for the 
clamped plates, the basis functions are given by 

^n(v) — cos(2 n - 1)tvt] - cos(2n + 1)7 rrj (10a) 

= cos2(n — 1)71-77 — cos2n7r?7 (10b) 

Algebraic equations governing buckling of long 
plates are obtained by substituting the series expan- 
sion for the buckling mode given by equation (8) 
into the second variation of the total potential en- 
ergy and then computing the integrals appearing in 
the second variation in closed form. The resulting 
equations constitute a generalized eigenvalue prob- 
lem that depends on the aspect ratio of the buckle 
pattern A/6 (see fig. 1) and the nondimensional pa- 
rameters defined by equations (1) through (4). The 
in plane loadings are expressed in terms of a load- 
ing parameter p that is increased monotonically un- 
til buckling occurs. The smallest eigenvalue of the 
problem p cr corresponds to buckling and is found by 
specifying a value of A/6 and solving the correspond- 
ing generalized eigenvalue problem for its smallest 
eigenvalue. This process is repeated for successive 
values of A/6 until the absolute smallest eigenvalue is 
found. The absolute smallest eigenvalue corresponds 
to p cr . 

Results obtained using the analysis described 
in the present paper were compared with results 


for isotropic and specially orthotropic plates pub- 
lished in the technical literature for several load- 
ing cases and boundary conditions. (See refs. 4 
through 10.) Results for anisotropic plates were also 
compared with results obtained using the computer 
code VIPASA. (See ref. 11.) In all comparisons, re- 
sults obtained from the analysis used in the present 
study were found to be within a few percent of the 
results published in references 4 through 10 and the 
VIPASA results. 

Laminate Construction and 
Nondimensional Parameters 

To characterize plate buckling behavior in terms 
of nondimensional parameters, it is important to un- 
derstand how laminate fiber orientation, stacking se- 
quence, and material properties affect their values. 
To accomplish this task and to determine the nu- 
merical range of the nondimensional parameters for 
typical balanced symmetric laminates, results are 
presented in this section for [(±#) m ] s angle-ply lam- 
inates, [(±45/0/90) m ] s and [(0/90/±45) m ] s quasi- 
isotropic laminates, and [(±45/02 ) m ] s and 

[(±45/902) m ] s laminates typically referred to as or- 
thotropic laminates. Each of these laminates pos- 
sesses specially orthotropic membrane deformation 
states prior to buckling. (That is, the A\§ and A26 
laminate constitutive terms are zero- valued.) 

Effects of Fiber Orientation and Stacking 

Sequence 

Values of the nondimensional parameters /?, 7, 6, 
and the parameter ( D 11 /D 22 )^ 4 arc shown in fig- 
ures 2 and 3 for [±1 9] s and [(±45) rn ] s angle-ply lam- 
inates, respectively. Similar results are presented 
in figure 4 for [(±45/0/90) J s and [(0/90/±45) m ] s 
quasi-isotropic laminates, and in figure 5 for 
[(±45/02) m ] s and [(±45/902) m ] s laminates. All re- 
sults presented in these figures are for laminates with 
plies made of a typical graphite-epoxy material hav- 
ing a longitudinal modulus E\ of 18.5 x 10 6 psi, a 
transverse modulus E 2 of 1.6 x 10 6 psi, an inplane 
shear modulus G12 of 0.832 x 10 6 psi, a major Pois- 
son’s ratio v \2 of 0.35, and a nominal ply thickness 
of 0.005 in. 

For the angle-ply laminates, the anisotropic pa- 
rameters 7 and 6 have maximum values at 9 = 39° 
and 51°, respectively (see fig, 2), and their values 
diminish monotonically as the number of plies in- 
creases (see fig. 3). Between the values of 39° and 
51°, values of 7 (and similarly <5) vary by only a 
couple of percent. The orthotropic parameters (3 
and {Du/ D22) 1 ^ have maximum values at 0 — 45° 
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and 0°, respectively, and their values remain con- 
stant as the number of plies increases. 

For the 0/45/90 family of quasi- isotropic lami- 
nates, the anisotropic parameters 7 and 8 shown in 
figure 4 asymptotically approach zero from above as 
the number of plies increases. The solid curves and 
dashed curves shown in figure 4 indicate results for 
the [(0/90/ ±45) m ] s and [(±45/0/90) TT? ] s laminates, 
respectively. The results in figure 4 indicate that the 
orthotropic parameters 0 and (D\\j for the 
[(±45/0/90) m ] s laminates asymptotically approach 
the value of 1 from above as the number of plies in- 
creases. The values of j3 and {D\i/ D 22 ) *^ 4 for the 
[(0/90/±45) m ] s laminates asymptotically approach 
the value of 1 from above and below, respectively, as 
the number of plies increases. A value of 1 for (3 and 
(T>ll/£>22) ly/4 anc * a va ^ uc °f 0 for 7 an( l ^ represent 
pure isotropy associated with a plate made from an 
essentially homogeneous material like wrought alu- 
minum plate. 

The results presented in figure 5 for the 
[(±45/02) m ] s and [(±45/902) m ] s laminates indicate 
that anisotropic parameters 7 and 6 also asymptoti- 
cally approach zero from above as the number of plies 
increases. The solid curves and dashed curves shown 
in figure 5 indicate results for the [( ±45/02 ) m ] s and 
[( ±45/902 ) m ] s laminates, respectively. The results in 
figure 5 indicate that both laminate families have the 
same value of the orthotropic parameter 0 for a given 
number of plies, and that (3 decreases monotonically 
to an asymptotic value as the number of plies in- 
creases. In contrast, the [(±45/0a) m ] s laminates have 

larger values of {D\\f 022 )^^ than the [(±45/902) m ] s 

laminates. Moreover, the values of (D\\/ D‘ 22 )^^ 
for the [(±45/90 2 ) m ] s and [(±45/02 ) m ] s laminates 
monotonically diminish and increase, respectively, as 
the number of plies increases. 

Effects of Varying Material Properties 

Results showing the effects of varying material 
properties on the nondimensional parameters for 
[±0] g angle-ply laminates are presented in figures 6 
through 9. In these figures, results are presented 
for laminates made of high-strength graphite-epoxy 
material, ultrahigh-modulus graphite-epoxy mate- 
rial, S-glass-epoxy material, Kevlar-epoxy material, 
boron-epoxy material, and boron-aluminum mate- 
rial. The moduli ratios E\jE<i and E2/G12 and the 
major Poisson’s ratio u \2 for these materials are given 
in table 1. 


Table 1. Stiffness Ratios and Major Poisson’s Ratio for 
Common Lamina Material Systems 


Material system 

i 

O 12 

v\i 

High-strength graphite-epoxy 

11.6 

1.9 

0.35 

Kevlar-epoxy 

13.8 

2.7 

0.34 

S-glass-epoxy 

4.4 

2,1 

0.25 

Ultrahigh-modulus graphite-epoxy 

50.0 

1.5 

0.26 

Boron-epoxy 

10.0 

4.3 

0.21 

Boron-aluminum 

1.6 

3.0 

0.23 


The results presented in figure 6 show that maxi- 
mum values of the orthotropic parameter 0 occur at 
9 = 45° for all the material properties considered. 
The largest values of 0 are, for most of the range of 
fiber orientation angles, exhibited by laminates made 
of ultrahigh-modulus graphite-epoxy material; these 
are followed by laminates made of Kevlar-epoxy ma- 
terial, boron-epoxy material, high-strength graphite- 
epoxy material, S-glass-epoxy material, and boron- 
aluminum material. The results presented in figure 7 
indicate that the largest value of ( D 11 /D 22 ) 1//4 is also 
exhibited by laminates made of ultrahigh-modulus 
graphite- epoxy material for 0° < 9 < 45°. The 
next highest values of (L>n/D22) 1 ^ 4 in the range of 
0 ° < 9 < 45° arc exhibited by laminates made of 
Kevlar-epoxy material, high-strength graphite-epoxy 
material, boron-epoxy material, S-glass-epoxy mate- 
rial, and boron- aluminum material. In the range of 
45° < 9 < 90°, the trend reverses. 

The results presented in figures 8 and 9 indi- 
cate that the largest values for the anisotropic pa- 
rameters 7 and 8 are also exhibited by laminates 
made of ultrahigh-modulus graphite-epoxy material 
followed by laminates made of Kevlar-epoxy mate- 
rial, boron-epoxy material, high-strength graphite- 
epoxy material, S-glass-epoxy material, and boron- 
aluminum material, i.e., the same trend as that 
exhibited by the orthotropic parameter 0. The 
maximum values of 7 and 8 for each material sys- 
tem range for 9 between 31° and 59°. For example, 
7 is at a maximum at 9 = 31° (7 = 0.13) and 8 is at a 
maximum at 9 = 59° (8 = 0.13) for laminates made 
of boron- aluminum material. For laminates made of 
ultrahigh-modulus graphite-epoxy material, 7 is at 
a maximum at 9 = 40° (7 = 0.68), and <5 is at a 
maximum at 9 — 50° (7 = 0.68). 

The results discussed in this section indicate a 
practical range of the numerical values of the non- 
dimensional parameters. In subsequent sections of 
this paper, generic buckling results are presented for 
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the ranges of 0.2 < (3 < 3.0, 0.5 < ( D 11 /D 22 < 
3.0, 0 < 7 < 0.8, and 0 < 6 < 0.8. Negative val- 
ues for 7 and 6 are possible for some laminates (e.g., 
[=p45] s laminates), but they are not considered in the 
present paper. An important physical consideration 
to keep in mind when using the nondimensional pa- 
rameters in an arbitrary manner to generate buck- 
ling design charts is that together they must sat- 
isfy the condition that the strain energy of the plate 
have a positive value. This condition is due to the 
thermodynamic restrictions that define the admissi- 
ble range of lamina material properties, and hence 
define the range of admissible plate bending stiff- 
nesses. This condition was enforced in the buckling 
analysis and accounts for the fact that some of the 
curves shown on the generic buckling design charts 
presented herein do not span the entire range of the 
nondimensional parameters specified herein. 

Buckling Results and Discussion 

Results are presented in this section for clamped 
and simply supported plates loaded by either ax- 
ial compression or shear. In addition, results 
are presented for plates subjected to a combined 
loading of axial tension or compression and shear, 
transverse tension or compression and shear, and 
transverse tension or compression and axial compres- 
sion. In each of the loading cases involving shear, a 
distinction between positive and negative shear 
loadings is made whenever anisotropy is present. A 
positive shear loading corresponds to a shear stress 
acting on the plate edge y = b (see fig. 1) in the 
positive x-coordinate direction. For each of the load- 
ing cases considered in the present paper, results are 
first presented for the familiar angle-ply laminates, 
and then generic results are presented in terms of the 
nondimensional parameters described herein. The 
results presented for the angle-ply laminates are cal- 
culated using lamina material properties of a typ- 
ical graphite-epoxy material. The properties used 
are the same ones that were used to calculate the 
nondimensional parameters shown in figure 2. 

Axial-Compression-Loaded Plates 

Nondimensional buckling coefficients K x as a 
function of the fiber orientation angle 9 arc presented 
in figure 10 for infinitely long, clamped and simply 
supported [±0] plates loaded in axial compression. 
The 4-ply [±0] s laminates are limiting cases, but they 
arc useful in the present study since they exhibit the 
highest degree of anisotropy of the angle-ply lami- 
nates. Similar results are presented in figure 11 for 
[(±45) m ] s plates for m — 1, 2, ..., 12. The solid 
curves shown in figures 10 and 11 correspond to re- 
sults for which anisotropy is included in the analysis. 


The dashed curves in the figures correspond to results 
for which the D\q and D2G anisotropic constitutive 
terms are neglected. The buckling coefficients shown 
in figure 10 were calculated for increments in 9 equal 
to 1°. Similar sets of results for square finite-length 
plates are presented in reference 1. 

The results presented in figure 10 indicate that ne- 
glecting plate anisotropy in the analysis always over- 
estimates the plate bending stiffness, and thus always 
yields nonconservativc estimates of the buckling re- 
sistance. The largest buckling coefficients for the or- 
thotropic plates occur at 9 — 45°, whereas the largest 
ones for the anisotropic plates occur at slightly less 
than 45°. The largest difference between the or- 
thotropic and anisotropic solutions for both clamped 
and simply supported plates is about 25 percent of 
the corresponding orthotropic solution. 

The results presented in figure 11 indicate that 
the anisotropic solutions for both clamped and sim- 
ply supported plates converge monotonically from 
below to the specially orthotropic solutions as the 
number of plies in the [(±45) m ] s laminates increases. 
The largest difference between the anisotropic and 
orthotropic solutions occurs for m = 1 and is approx- 
imately 25 percent of the corresponding orthotropic 
solution. The difference is approximately 5 percent 
for m = 2 and less than 1 percent for m = 6. 

Generic buckling results for clamped compression- 
loaded plates are presented in figures 12 through 14 
for a wide range of values of the four nondimensional 
parameters. The results presented in figure 12 show 
buckling coefficients of clamped plates as a function 
of the parameter (Dn/Z^) 1 ^ 4 for P = 3. Buckling 
curves are presented in this figure for discrete values 
of 7 = 6 ranging from 0 to 0.8. The [±45] s lami- 
nates made of a typical graphite-epoxy material and 
an ultrahigh-modulus graphite- epoxy material have 
values of 7 and 6 equal to 0.52 and 0.68, respec- 
tively. Thus, a value . of 0.8 corresponds to an ex- 
tremely anisotropic plate. The straight lines shown 
in figure 12 are parallel to the abscissa axis and thus 
indicate that the buckling coefficients are indepen- 
dent of the parameter {D \\/ for the entire 
range of anisotropy shown in the figure. 

Similar plots were made for clamped and sim- 
ply supported plates in which the orthotropic pa- 
rameter P was varied and unequal values of 7 and 
<5 were used. The additional results also show the 
buckling coefficients to be independent of the pa- 
rameter (D\i/ T>22) 1//4 - This finding is consistent 
with the observation that plots of buckling coeffi- 
cient versus plate aspect ratio (length over width) 
for finite-length plates attenuate to constant values 
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as the plate aspect ratio becomes large. More specif- 
ically, in a finite-length plate the bending stiffness 
ratio (On/^22) 1 ^ 4 influences the size and number of 
buckles that occur along the plate length. The num- 
ber of buckles that occur along the plate length, in 
turn, directly affects the value of the buckling coef- 
ficient for a given plate (as indicated by the festoon 
nature of the buckling curves for finite-length plates). 
For an infinitely long plate, the bending stiffness ra- 
tio (D11/D22) 1 / 4 still influences the size of the buck- 
les with respect to the basic periodic unit of an in- 
finitely long plate, but the number of buckles that 
occur along the plate length becomes meaningless. 

The independence of K x with respect to 
(£>u/£>22) 1/4 represents an important simplification 
in that the buckling coefficients of long plates can 
be represented by a single orthotropic parameter, 
namely, (3. However, the parameter (^n/^22) ly/4 is 
important since it does affect the aspect ratio of the 
buckle pattern, and it is useful in determining when 
classical plate theory becomes insufficient and trans- 
verse shear deformation must be included in buckling 
analyses. 

Buckling coefficients arc presented in figure 13 for 
clamped compression- loaded plates as a function of 
the orthotropic parameter (3 for values of 7 = 8 rang- 
ing from 0 to 0.8. Using the result that the buck- 
ling coefficients are independent of (Dn/Z^) 1 ^ 4 , re- 
sults corresponding to particular laminates made of 
a typical graphite-epoxy material and an ultrahigh- 
modulus graphite- epoxy material are also shown in 
figure 13 and are indicated by symbols. The symbols 
that represent several laminates indicate that all the 
laminates have the same buckling coefficient. The re- 
sults presented in this figure indicate that anisotropy 
can significantly reduce buckling resistance and that 
the effect of anisotropy tends to diminish somewhat 
as (3 increases (i.c., the curves become spaced closer 
together as (3 increases). Similar results were ob- 
tained for simply supported plates. These results in- 
dicate the same trends and indicate that the effects 
of varying j3 are slightly more pronounced for the 
clamped plates. (That is, the linear portions of the 
curves for the clamped plates generally had a higher 
slope than the corresponding curves for the simply 
supported plates.) 

Results are presented in figure 14 for clamped 
compression-loaded plates in which the anisotropic 
parameters 7 and 8 are varied independently. Curves 
are also presented in this figure for values of 7 and 
8 ranging from 0 to 0.8. Moreover, two sets of 
curves are shown in figure 14 corresponding to values 
of the orthotropic parameter (3 equal to 1 and 3. 


The results presented in this figure indicate that 
the buckling coefficients are more sensitive to the 
parameter 8 than to 7. This observation can be 
seen by comparing the differences in the values of the 
buckling coefficients in which 7 is held constant and 

8 is varied, and vice versa. The results also support 
the previous statement that the effect of anisotropy 
tends to diminish somewhat as (3 increases. A set 
of results similar to those presented in figure 14 
were obtained for simply supported plates. These 
results indicate the same trends and indicate that 
the effects of varying /?, 7, and 8 arc also slightly 
more pronounced for the clamped plates. 

Shear-Loaded Plates 

Shear buckling coefficients K s arc presented in 
figure 15 as a function of the fiber orientation angle 

9 for infinitely long, clamped and simply supported 
[±0] s plates. Similar results are presented in figure 16 
for [(±45) m ] s plates for m = 1 , 2, . . . , 12. The solid 
and dotted curves shown in figures 15 and 16 corre- 
spond to results for which anisotropy is included in 
the analysis. The solid curves correspond to a posi- 
tive shear loading, and the dotted curves correspond 
to a negative shear loading. The dashed curves in the 
figures correspond to results for which anisotropy is 
neglected, and no distinction between positive and 
negative shear is necessary. The buckling coefficients 
shown in figure 15 were calculated for increments in 
9 equal to 1°. 

The results presented in figure 15 indicate that ne- 
glecting plate anisotropy in the analysis yields non- 
conservative estimates of the buckling resistance of 
plates loaded in positive shear and conservative esti- 
mates for the plates loaded in negative shear. More- 
over, the plates loaded in negative shear exhibit the 
most buckling resistance. The largest buckling coeffi- 
cients occur at 9 = 45° for the orthotropic plates and 
at slightly more than 45° for the anisotropic plates 
loaded in negative shear. The largest buckling co- 
efficients occur at 9 = 0° and 90° for the clamped 
plates loaded in positive shear and at 20° for the cor- 
responding simply supported plates. The largest dif- 
ference between the orthotropic and anisotropic solu- 
tions for both clamped and simply supported plates 
is approximately 48 percent of the corresponding or- 
thotropic solution. This difference is nearly twice 
that of the corresponding compression-loaded plates. 

The results presented in figure 16 indicate that 
the anisotropic solutions for both clamped and sim- 
ply supported plates converge monotonically to the 
specially orthotropic solutions as the number of plies 
increases in the [(±45) m ] s laminates. The plates 
loaded by negative shear converge from above, and 
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the plates loaded by positive shear converge from be- 
low. The largest difference between the anisotropic 
and orthotropic solutions occurs for m — 1 and is 
approximately 48 percent of the corresponding or- 
thotropic solution. The difference is approximately 
7 percent for m = 6 and approximately 3.5 percent 
for m = 12. Comparing the results presented in fig- 
ures 11 and 16 indicates that the effects of anisotropy 
are much more pronounced in shear-loaded plates 
than in compression- loaded plates. 

Generic buckling results for clamped and simply 
supported shear-loaded plates arc presented in fig- 
ures 17 through 19 for a wide range of values of the 
nondimensional parameters. The results presented 
in figure 17 show buckling coefficients of clamped 
plates as a function of the parameter {D\\/ 7^22 ) 
for P = 3. Buckling curves are presented in this fig- 
ure for values of 7 = 6 ranging from 0 to 0.8. Like 
the compression-loaded plates, the results presented 
in figure 17 indicate that the buckling coefficients 
are independent of the parameter (D\\/ for 
the entire range of anisotropy. Similar results were 
obtained for clamped and simply supported plates 
in which the orthotropic parameter (3 was varied. 
The additional results also show that the buck- 
ling coefficients are independent of the parameter 
(Dii/T? 22) 1 ^' The parameter (Dn/D22) 1/A4 was 
found, however, to affect the aspect ratio of the 
buckle pattern. 

Buckling coefficients are presented in figure 18 for 
clamped shear-loaded plates as a function of the or- 
thotropic parameter /?. Curves are also presented 
in this figure for values of 7 = 6 ranging from 0 
to 0.8. Since the buckling coefficients are inde- 
pendent of (Z)u/D22) 1 ^ 4 , the results corresponding 
to particular laminates that are made of a typical 
graphite-epoxy material and an ultrahigh-modulus 
graphite-epoxy material are indicated in the figure by 
symbols. Like the compression- loaded plates, the re- 
sults presented in this figure indicate that anisotropy 
can significantly reduce buckling resistance of shear- 
loaded plates and that the effects of anisotropy tend 
to diminish somewhat as (3 increases. Similar results 
were obtained for simply supported plates. These re- 
sults also indicate the same trends and demonstrate 
that the effects of varying (3 are slightly more pro- 
nounced for the clamped plates. 

Results are presented in figure 19 for clamped 
shear-loaded plates as a function of the anisotropic 
parameters 7 and 6 with values ranging from 0 
to 0.8 for /3 == 3. The results presented in figure 19 
indicate the same trend shown in figure 14 for the 
compression-loaded plates; i.e., the buckling coeffi- 


cients are more sensitive to 6 than to 7. Similar re- 
sults were obtained for (3 = 1. These results also 
indicate that the effects of anisotropy tend to dimin- 
ish somewhat as (3 increases. Results similar to those 
presented in figure 19 were obtained for simply sup- 
ported plates with (3=1 and 3. These results indi- 
cate the same trends and once again reveal that the 
effects of varying /?, 7, and 6 are slightly more pro- 
nounced for the clamped plates than for the simply 
supported plates. 

Plates Subjected to Combined Loadings 

Results are presented in this section for plates 
subjected to combined loadings. First, results are 
presented for plates subjected to axial tension or 
compression loadings combined with shear loading. 
Next, results are presented for plates subjected to 
transverse tension or compression loadings combined 
with shear loading. Last, results are presented for 
plates subjected to axial compression combined with 
either transverse tension or compression loadings. 

Axial tension or compression loadings 
combined with shear loading. Buckling inter- 
action curves are presented in figure 20 for infinitely 
long clamped [±9] s plates subjected to axial tension 
or compression and shear loadings for several fiber 
orientation angles 6. The dashed curves shown in 
figure 20 correspond to results for which anisotropy 
is included in the analysis, and the solid curves 
correspond to results for which the D\§ and D26 
anisotropic constitutive terms are neglected. Neg- 
ative values of K x correspond to axial tension load- 
ings, and negative values of K s correspond to neg- 
ative shear loadings. Points of intersection of the 
interaction curves shown in figure 20 with a straight 
line emanating from the origin of the plot constitute 
constant values of ( N xy /N x ) (Du /D 22 ) and are 
depicted in the figure for convenience. 

The results presented in figure 20 indicate that 
the effects of anisotropy manifest themselves as a 
phase shift of self-similar curves in the AT- Airplane. 
The self-similar nature of the curves suggests that it 
may be possible to obtain simple buckling interaction 
formulas like those existing for isotropic plates. (See 
ref. 5.) For plates loaded in positive shear, neglecting 
plate anisotropy in an analysis yields buckling coef- 
ficients substantially larger than those correspond- 
ing to the anisotropic solution. For plates loaded in 
negative shear, however, this trend is reversed. Sim- 
ilar results were obtained for corresponding simply 
supported plates, and these results exhibit exactly 
the same trends. In addition, both the results in 
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figure 20 and the corresponding results for simply 
supported plates suggest that the buckling resistance 
of a compression-loaded plate can be improved by ap- 
plying a negative shear loading. For the plates with 
9 — 45 0 5 the improvement is approximately 20 per- 
cent and 25 percent of the K x values corresponding to 
pure axial compression loading for the clamped and 
simply supported plates, respectively. The results 
presented in figure 20 also indicate that laminates 
with 6 — 45° yield the most buckling-resistant inter- 
action curves for both cases in which anisotropy is 
included and neglected. This trend is also exhibited 
by the simply supported plates. 

Generic buckling results for clamped and sim- 
ply supported plates subjected to combined axial 
and shear loadings are also presented in figures 21 
through 23 for a wide range of values of the non- 
dimensional parameters. Additional results were ob- 
tained for both boundary condition cases and several 
values of /?, 7, and 6 . These additional results in- 
dicate that the buckling interaction curves, such as 
those presented in figures 21 through 23, are inde- 
pendent of the parameter (D\ 1/D22) 

Buckling interaction curves are presented in fig- 
ure 21 for clamped and simply supported plates 
and for several values of the orthotropic parame- 
ter p. The curves appearing in this figure arc for 
values of 7 = 6 — 0. By using the simplification 
that the buckling interaction curves are independent 
of (Z?h/Z? 22 ) 1 ^j results corresponding to particular 
laminates can be obtained once the value of P for the 
specific laminate is known, provided the correspond- 
ing nondimensional anisotropic parameters have a 
value of zero. Good estimates of buckling resistance 
can be otained for laminates with small values of the 
nondimensional anisotropic parameters compared to 
unity. For example, the curves with P = 1 correspond 
to solutions for isotropic plates and give good ap- 
proximations to the results for [(±45/0/90) m ] g quasi- 
isotropic laminates with 32 or more plies made of a 
typical graphite-epoxy material. An important result 
shown in figure 21 is that substantial gains in buck- 
ling resistance can be obtained by tailoring laminate 
construction to increase the parameter p. 

Results are presented in figure 22 for clamped 
and simply supported plates with values of 7 = <5 
ranging from 0 to 0.8 and for /? — 3. These results 
also show that anisotropy is generally manifested as 
a phase shift of the interaction curves in the K X K S - 
plane. The results in figure 22 support the state- 
ment given previously that the buckling resistance 
of a compression-loaded plate can be improved by 
applying a negative shear loading (as was seen for 


the angle- ply laminates). For extremely anisotropic 
plates with 7 = 6 = 0.8, the improvement indi- 
cated by these results is approximately 90 percent 
and 85 percent of the K x values corresponding to 
pure compression loading for the clamped and sim- 
ply supported plates, respectively. The results pre- 
sented in figure 22 also indicate that plates with 
7 = 6 = 0 are the most buckling resistant of the 
plates loaded by positive shear forces, whereas the 
plates with 7 = 6 = 0.8 are the most buckling re- 
sistant of the plates loaded by negative shear forces. 
This trend is also exhibited by the simply supported 
plates. 

Results are presented in figure 23 for clamped 
plates with P — 1 as a function of the anisotropic 
parameters 7 and <5. Curves are shown in this 
figure for values of 7 = 0 and 0.4 and for values 
of 6 ranging from 0 to 0.6. The results presented 
in figure 23 indicate the same trend as that shown 
in figure 14 for the axial compression-loaded plates 
and in figure 19 for the shear-loaded plates; i.c., 
the buckling coefficients are more sensitive to the 
anisotropic parameter 6 than to 7. 

Transverse tension or compression load - 
ings combined with shear loading . Buckling 
interaction curves are presented in figure 24 for in- 
finitely long clamped [±0] s plates subjected to trans- 
verse tension or compression and shear loadings for 
several fiber orientation angles 9. The dashed curves 
shown in figure 24 correspond to results for which 
anisotropy is included, and the solid curves corre- 
spond to solutions for which it is neglected. Neg- 
ative values of K y correspond to transverse tension 
loadings, and negative values of K s correspond to 
negative shear loadings. Points of intersection of the 
interaction curves shown in figure 24 with a straight 
line emanating from the origin of the plot constitute 
constant values of (N y /N xy ) (D\\f D22) 1 ^ • The hor- 
izontal straight-line portions of the buckling inter- 
action curves indicate values of the buckling coeffi- 
cients for which a plate buckles into a wide column 
mode (Ky = 4 for clamped plates). For this mode, 
the half-wavelength A shown in figure 1 approaches 
infinity in the analysis. 

The results presented in figure 24 indicate that 
the effects of anisotropy also manifest themselves as 
a phase shift of self-similar curves in the A^ Ay-plane 
for these loading conditions. For plates loaded in 
positive shear, neglecting plate anisotropy in an anal- 
ysis yields buckling coefficients substantially larger 
than those corresponding to the anisotropic solu- 
tion. In addition, the value of K s at which the 
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plate buckles into a wide column mode is significantly 
overestimated. For plates loaded in negative shear, 
this trend is reversed. The coalescence of all the 
curves shown in figure 24 into the same horizontal- 
line segment at K y = 4 indicates that anisotropy 
is unimportant for plates that buckle into a wide 
column mode; i.e., the transverse loading N y forces 
the plate to assume a deformation state that is not 
skewed (or it is skewed at x = ±oc). Similar results 
were obtained for corresponding simply supported 
plates, and these results exhibit exactly the same 
trends. The simply supported plates buckle into a 
wide column mode at K y = 1. 

Generic buckling results for clamped and simply 
supported plates subjected to combined axial and 
shear loadings are presented in figures 25 through 27. 
Additional results were obtained for both boundary 
condition cases and several values of /?, 7, and 5. 
These additional results indicate that the buckling 
interaction curves for this combination of loadings 
are also independent of the parameter {D\i/D22) 1 ^- 

Buckling interaction curves are presented in fig- 
ure 25 for clamped and simply supported plates 
with several values of the orthotropic parameter (3 
and with 7 = 5 = 0. Once again, the fact that 
the buckling interaction curves are independent of 
(Dii/T> 22) 1//1 permits results for particular lami- 
nates with zero-valued (or small compared to unity) 
anisotropic parameters to be obtained directly from 
figure 25 once the value of (3 for the specific lami- 
nate is known. The results shown in figure 25 in- 
dicate that substantial gains in buckling resistance 
can be obtained by tailoring laminate construction 
to increase the parameter /?. 

Generic buckling results arc presented in fig- 
ure 26 for clamped and simply supported plates with 
values of 7 = 8 ranging from 0 to 0.8 and for a value 
of P = 3. These results also confirm the results pre- 
sented in figure 24 that show that the anisotropy is 
generally manifested as a phase shift of the inter- 
action curves in the K s K y ~ plane, and that anisotropy 
does not affect the wide column buckling mode. The 
results presented in figure 26 also indicate that plates 
with 7 = 5 = 0 are the most buckling resistant of the 
plates loaded by positive shear forces, whereas the 
plates with 7 = 5 = 0.8 are the most buckling resis- 
tant of the plates loaded by negative shear forces. 

Results are presented in figure 27 for clamped 
plates with /? = 3 in which the anisotropic param- 
eters 7 and 5 are varied independently. Buckling 
interaction curves arc shown in this figure for 7 = 0 
and 0.6 and for values of 5 ranging from 0 to 0.6. 
The results presented in figure 27 indicate that the 


plates loaded by transverse tension or compression 
and shear are also more sensitive to the anisotropic 
parameter 5 than to 7. 

Axial compression combined with either 
transverse tension loading or compression 
loading . Buckling interaction curves are presented 
in figure 28 for infinitely long, clamped and simply 
supported [±0] s plates subjected to transverse ten- 
sion or compression and axial compression loadings 
for several fiber orientation angles 6. The dashed 
curves shown in figure 28 correspond to results for 
which anisotropy is included in the analysis, and the 
solid curves correspond to results for which it is ne- 
glected. In this figure the points of intersection of 
the interaction curves with a straight line emanating 
from the origin of the plot constitute constant values 
of {Ny/Nx) (D\\fD 22 ) 1 ^ 2 - The horizontal straight- 
line portions of the buckling interaction curves in- 
dicate values of the buckling coefficients at which a 
plate buckles into a wide column mode. 

The results presented in figure 28 indicate that 
the effects of anisotropy for plates subjected to these 
loading conditions also generally manifest themselves 
as a phase shift of self-similar curves in the K x K y - 
plane. Neglecting plate anisotropy in an analysis 
always yields buckling coefficients larger than those 
corresponding to the anisotropic solution. Similar to 
the plates loaded by shear and transverse compres- 
sion, the value of K x at which the plate buckles into 
a wide column mode is substantially overestimated 
when anisotropy is neglected. 

Generic buckling results are presented in fig- 
ures 29 through 31 for clamped and simply sup- 
ported plates subjected to combined axial and trans- 
verse loadings. Additional results were obtained for 
both boundary condition cases and several values of 
,3, 7, and 5. These additional results indicate that 
the buckling interaction curves for this combination 
of loadings are also independent of the parameter 
(Dn/^22) 1 ^. 

Buckling interaction curves are presented in fig- 
ure 29 for clamped and simply supported plates with 
several values of the orthotropic parameter /? and 
with 7 = 5 = 0. Since the buckling interaction curves 
are independent of {D n / D 2 2 ) ] ^ , results for particu- 
lar laminates with zero- valued (or small compared to 
unity) nondimensional anisotropic parameters can be 
obtained directly from figure 29 once the value of j3 
for the specific laminate is known. The results shown 
in figure 29 indicate that substantial gains in buck- 
ling resistance can be obtained when (3 increases. 
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Generic buckling results are presented in figure 30 
for clamped plates with values of 7 = 8 ranging from 
0 to 0.6. Two sets of curves are shown in the figure. 
The solid and dashed curves correspond to values of 
(3 = 3 and 1, respectively. These results also show a 
phase shift due to anisotropy, and they show that the 
effects of anisotropy are more pronounced for plates 
with (3 = 1 than for plates with (3 = 3. The results 
presented in figure 30 also indicate that plates with 
7 — £ — 0 are the most buckling resistant of the 
plates when nonzero axial loading is present. 

Results arc presented in figure 31 for clamped 
plates with (3 — 3 as a function of the anisotropic 
parameters 7 and 8. Curves are shown in this figure 
for 7 = 0 and 0.6 and for values of 8 ranging from 0 to 
0.6. The results presented in figure 31 indicate that 
the plates loaded by transverse tension or compres- 
sion and axial compression are also generally more 
sensitive to the parameter 8 than to 7. 

Concluding Remarks 

A parametric study has been presented of the 
buckling behavior of infinitely long, symmetrically 
laminated anisotropic plates subjected to combined 
loadings. Loading conditions consisting of axial ten- 
sion and compression, transverse tension and com- 
pression, and shear were investigated for clamped 
and simply supported plates. Results are presented 
that were obtained using a special-purpose analysis 
that was derived in terms of useful nondimensional 
parameters. The analysis was found to be well-suited 
for parametric studies. Buckling results are also 
presented for some common laminate constructions, 
and generic buckling design charts have been pre- 
sented for a wide range of parameters. The generic 
design charts are presented in terms of useful non- 
dimensional parameters, and the dependence of the 
nondimensional parameters on laminate fiber orien- 
tation, stacking sequence, and material properties is 
discussed. 

The results presented in the present, paper show 
that nondimensional parameters can be very use- 
ful in presenting results in a concise manner for a 
wide range of loading conditions, boundary condi- 
tions, and laminate constructions. Results are pre- 
sented that show that, the effects of anisotropy are 
much more pronounced in shear-loaded plates than in 
compression-loaded plates. In addition, the effects of 
anisotropy on plates subjected to combined loadings 
are shown to be generally manifested as phase shifts 
of self-similar buckling interaction curves. These re- 
sults indicate that the buckling resistance of long 
compression-loaded, highly anisotropic plates can be 
improved significantly by applying a shear loading 


with a specific orientation. Moreover, it is shown 
that anisotropy reduces the buckling resistance of bi- 
axially compressed plates. 

A substantial number of generic buckling re- 
sults are presented in the present paper for a wide 
range of values of the nondimensional parameters. 
In all cases considered in the study, the buck- 
ling coefficients of infinitely long plates arc shown 
to be independent of the bending stiffness ratio 
(Dn/I>22) 1/ ^ 4 * Also, large increases in buckling 
resistance are shown to be obtained by tailoring 
the laminate construction to increase the parameter 
p = (D12 + 2D$(>) / (D11D22) 1 / 2 , and the importance 
of anisotropy generally diminishes as (3 increases. 
Results are also presented that show thc buckling 
coefficients to be generally more sensitive to the 
anisotropic parameter 8 = D26/ (R11 ^22) ^ than to 
the anisotropic parameter 7 = 1 ^22) ly/4 for 

the entire range of loadings and boundary conditions 
considered. 


NASA Langley Research Center 
Hampton, VA 23665-5225 
March 25, 1992 
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Figure 1. Geometry of buckled plate subjected to biaxial compression and shear loadings. 
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Figure 2. Nondimensional parameters for [±0] s graphite-epoxy laminates. 
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Figure 3. Nondimcnsional parameters for [(±45) m ] s graphite-epoxy laminates. 
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Figure 4. Nondimensional parameters for quasi-isotropic graphite-epoxy laminates. 
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Figure 5. Nondimensional parameters for [(±45/0 2 ) m ] s and [(±45/90 2 )J s graphite-epoxy laminates. 
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Figure 6. Effect of lamina material properties on nondimensional orthotropic parameter ft for [±0] s laminates. 
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Figure 7 . Effect of lamina material properties on 


parameter {D\y/ D 2 2 ) 1 / 4 for [± 0 ] s laminates. 


♦ y Material 

High-strength graphite-epoxy 


q 1 Kevlar-epoxy 

f J — ► x S-glass-epoxy 



Figure 8. Effect of lamina material properties on nondimensional anisotropic parameter 7 for [±0] g laminates. 
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Figure 9. Effect of lamina material properties on nondimensional anisotropic parameter 5 for [±#] s laminates. 
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Figure 10. Buckling coefficients for axial-compression-loaded [±#] s graphite-epoxy laminates. 
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Figure 11. Buckling coefficients for axial-compression-loaded [(±45) m ] s graphite-epoxy laminates. 
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Figure 12. Effects of (Du/D22) 1 ^ 4 and anisotropic parameters 7 and 6 on buckling coefficients for axial- 
compression- loaded clamped plates for P = 3. 
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Figure 13. Effects of orthotropic parameter j3 and anisotropic parameters 7 and 8 on buckling coefficients for 
axial-compression-loaded clamped plates. All laminates shown arc made of typical graphite-epoxy material 
unless otherwise noted. 
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Figure 14. Effects of anisotropic parameters 7 and 8 on buckling coefficients for axial- compression-loaded 
clamped plates. 
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Figure 15. Buckling coefficients for shear-loaded [±0] s graphite-epoxy laminates. 
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Figure 17. Effects of (D\\/ i? 22)^ 4 and anisotropic parameters 7 and 6 on buckling coefficients for shear-loaded 
clamped plates for /3 = 3. 
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Figure 18. Effects of orthotropic parameter /? and anisotropic parameters 7 and <5 on buckling coefficients 
for shear-loaded clamped plates. All laminates shown are made of typical graphite-epoxy material unless 
otherwise noted. 
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Figure 19. Effects of anisotropic parameters 7 and 6 on buckling coefficients for shear-loaded clamped plates 
for (3 = 3 . 



Figure 20. Buckling interaction curves for clamped [±0] g graphite-epoxy laminates subjected to axial tension 
or compression and shear. 



Figure 21. Effect of orthotropic parameter (3 on buckling interaction curves for specially orthotropic laminates 
(^ = <5 = 0 ) subjected to axial tension or compression and shear. 



Figure 22 . Effect of anisotropic parameters 7 and 6 on buckling interaction curves for laminates subjected to 
axial tension or compression and shear for /? = 3. 
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Figure 23. Effect of anisotropic parameters 7 and <5 on buckling interaction curves for clamped laminates 
subjected to axial tension or compression and shear for j3 = 1 . 



Figure 24. Buckling interaction curves for clamped [±0] s graphite-epoxy laminates subjected to transverse 
tension or compression and shear. 
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Figure 26. Effect of anisotropic parameters 7 and S on buckling interaction curves for laminates subjected to 
transverse tension or compression and shear for 0 = 3. 
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Figure 27. Effect of anisotropic parameters 'y and 5 on buckling interaction curves for clamped laminates 
subjected to transverse tension or compression and shear for /? = 3. 
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Figure 28. Buckling interaction curves for [±#] s graphite-epoxy laminates subjected to transverse tension or 
compression and axial compression. 
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Figure 29. Effect of orthotropic parameter /? on buckling interaction curves for specially orthotropic laminates 
(7 = $ — 0 ) subjected to transverse tension or compression and axial compression. 



Figure 30. Effect of anisotropic parameters 7 and 6 on buckling interaction curves for clamped laminates 
subjected to transverse tension or compression and axial compression. 
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Figure 31. Effect of anisotropic parameters 7 and 6 on buckling interaction curves for clamped laminates 
subjected to transverse tension or compression and axial compression for (3 — 3 . 
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